Electronic Structure of CeRhlns: dHvA and Energy Band Calculations. 
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The de Haas - van Alphen effect and energy band calculations are used to study angular dependent 
extremal areas and effective masses of the Fermi surface of the highly correlated antiferromagnetic 
material CeRhlns. The agreement between experiment and theory is reasonable for the areas mea- 
sured with the field applied along the (100) axis of the tetragonal structure, but disagree in size 
for the areas observed for the field applied along the (001) axis where the antiferromagnetic spin 
alignment is occurring. Detailed comparisons between experiment and theory are given. 
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I. INTRODUCTION 

The compounds CeMIus (M=Co, Ir, Rh) are a ne 
reported family of heavy fermion superconductorsEI" 
These materials crystallize in the tetragonal HoCoGas 
structure and are built of alternating stacks of Celna 
and MIn2. CeCoIns and Celrlns have superconducting 
transition temperatures of 2.3 K and 0.4 K, respectively, 
whereas CeRhlns orders antiferromagnetically at 3.8 K 
at ambient pressure. Applied pressures of order 16 kbar 
can induce an apparently first-order transition from the 
magnetically ordered state to a superconducting one with 
Tc= 2.1 K. The particular attraction of these materials 
is that not only do they represent a substantial increase 
in the number of known heavy fermion superconductors 
but also they appear to be quasi-2D variants of Celns , an 
ambient-pressure antiferromagnet in which stuperconduc- 
tivity can be induced at 25 kbar and 100 mK.B If one can 
demonstrate that the reduced dimensionality is respon- 
sible for the factor of ten increase in superconducting Tj, 
as well as the anomalous evolution from antiferromagnet 
to superconductor with pressure, their impact will far 
exceed that of just another family of heavy fermion su- 
perconductors. However, much work remains to establish 
this hypothesis. 

Because it displays both 'high-temperature' supercon- 
ductivity as well as unconventional magnetic behavior, 
CeRhlns seems an ideal candidate for electronic and mag- 
netic spectroscopic probes to establish the degree of re- 
duced dimensionality and its impact on fluctuation spec- 
tra. Some progress has been Wgde in establishing the 
magnetic structure of CeRhlnsBQ. Nuclear quadrupole 
resonance and neutron diffraction measurements have 
shown that the magnetically ordered state of CeRhlns is 
a spiral spin structure in which the spins in a given Celna 
plane are ordered antiferromagnetically, and the direction 
of the ordered moments in the plane spirals along the c- 
axis, rotating 107 degrees per layer. Although complete 



inelastic neutron scattering measurements of the dynam- 
ical susceptibility remain in progress, already NQR has 
established that the temperature evolution of the sub- 
lattice magnetization below T^r is much more rapid than 
would be expected from mean field theory and is a sign 
of reduced dimensionality. 

To date, studies of the electronic structure c£ CeRhlns 
have been far less complete. Cornelius et aln have re- 
ported limited de Haas-van Alphen (dHvA) measure- 
ments along the principal crystallographic axes using 
pulsed fields. These measurements reveal an anomalous 
temperature-dependence of the dHvA amplitude as well 
as large directional anisotropics in the effective masses of 
carriers that have been observed. 

Here, we report a more comprehensive dHvA study 
throughout the entire Brillouin zone of CeRhlns 
and compare our results directly with local-density- 
approximation band structure calculations. -Recently, a 
similar study of Celrlns has been reported.EI Taken to- 
gether, our results and those of Reference ^ allow di- 
rect evaluation of the changes in electronic properties of 
CeMIns as a function of transition metal ion M. 

In what follows, we report the details of our dHvA 
measurements and band structure calculations, discuss 
the extent to which they are self-consistent, and compare 
and contrast the electronic structures of CeRhlns and 
Celrlns. 



II. MEASUREMENTS 

All of the measurements reported here were made at 
the National High Magnetic Field Laboratory, Tallahas- 
see, FL using cantilever magnetometry at temperatures 
between 20 and 500 mK in applied fields ranging from 5 
to 18 T. Complete field rotations in the [100] and [001] 
planes of the tetragonal structure are reported. The sam- 
ple was grown from an In flux and etched in a 25% HCl in 
H2O solution down to a small plate that was mounted on 
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the cantilever with Apiezon N vacuum grease or GE Var- 
nish. This technique is most sensitive to low frequency 
oscillations not observed typically with the magnetic field 
modulation technique. 

The dHvA effect, a measurement of the oscillatory part 
of the magnetization, is a method for determining Fermi 
surface properties. The oscillatory magnetization M is 
given by the Lifshitz-Kosevitch (LK) equation (see Ref. 

for the mathematical details): 



M = -2.602 X 10" 
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where a — 1.47(m/mo) x 10^ G/K, A" is the second 
derivative of the area of the Fermi surface (FS) cross- 
section that is perpendicular to the applied field, G is the 
spin reduction factor, p is the harmonic number, and x 
is the Dingle temperature. The frequencies of the dHvA 
oscillations are proportional to extremal areas of the FS, 
and Fermi-liquid theory works well for heavy- fecmion ma- 
terials as has been shown in previous studies.Ell 

Since the magnetization of the sample was measured 
with a torque cantilever, it is necessary to understand 
how torque is related to the LK equation. dHvA oscil- 
lations in the torque arise from anisotropy in the Fermi 
surface, such that 
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where F is the dHvA frequency, 6 is the angle of the ap- 
plied field, M is the LK expression above, and V is the 
volume of the sample. This means that a roughly spher- 
ical FS will have a smaller torque signal than a highly 
elliptical FS. 

Two versions of a Fourier transform algorithm were 
used in the analysis: preliminary analysis was done with 
a fast Fourier transform (FFT) and in some cases a dis- 
crete Fourier transform (DFT) was used to increase the 
frequency resolution. Fourier spectra of the oscillations 
obtained for the field parallel to [100] and [001] are shown 
in Figure ^ and Figure |^. As can be seen there are sev- 
eral fundamental frequencies and several combinations 
(especially with the F7; e.g. F7 - Fg = 1302 T). Similar 
combination frequencies were found in all field directions; 
their origins will be discussed below. 

The fundamental frequencies are plotted as a func- 
tion of angle in Figure ||. As can be seen the mea- 
sured frequencies span the range from about 100 T over 
wide angular ranges to 5 kT for the field applied near 
[001]. The one defect of the cantilever technique was al- 
luded to in the previous paragraph - small ^ produces a 
small signal amplitude. The measured dHvA oscillations 
were small for [001] and [100] orientations, and the sig- 
nal amplitude was maximized at approximately a forty- 
five degree rotation from each of these principal axes. 



One might propose that field modulation measurements 
would be more appropriate in this case, but given the low 
(less than 1 kT) frequencies measured here, the required 
modulation fields are technically impractical in the high 
field ranges used. 



III. DISCUSSION 

The total energy and band structure of CeRhlns 
were calculated with a full-potential electronic structup^ 
method that uses linear mufHn-tin orbitals as bases.llHl 
Calculations were performed to obtain theoretical struc- 
tural parameters, one-electron bands and spectral densi- 
ties, and the structure of the Fermi surface. All calcula- 
tions used the GGAt^ to treat exchange and correlation 
and included the spin-orbit interaction in the variational 
basis. In the range of parameters searched, no magnetic 
instabilities were found. The calculations reported here 
are for paramagnetic CeRhlns in the simple tetragonal 
structure (P4/mmm). 

Structural parameters were obtained by simultaneous 
variation of volume, c/a, and the single internal struc- 
tural parameter to minimize the total energy. Substan- 
tially converged results were obtained with a set of 45 
irreducible (486 total) Brillouin zone points. The cal- 
culated volume, c/a ratio, and internal parameter are 
given in Table ||. Conventional electronic structure the- 
ory (with itinerant Ce 4/ electrons) applied to a simple 
paramagnetic cell seems to give a satisfactory descrip- 
tion of the structural properties of CeRhlns . Also given 
in Table | are the bulk modulus calculated at both the 
experimental and theoretical volumes and the Density of 
States (DOS) at the Fermi energy (Ep). calculated at 
the experimental volume. 

The band structure, calculated at the experimental 
volume, along high symmetry lines in the Brillouin zone 
is shown in Figure 0. Bands of predominantly Ce 4/ 
character dominate the band structure at and just above 
Fip. The set of bands along X T ^ Z ^ R are indica- 
tive of the nesting found in the Fermi surface discussed 
below. The electron DOS, shown in Figure |^, shows the 
predominance of Ce 4/ character at the Fermi energy. 
The DOS at Fip, while fairly large, does not indicate a 
dramatic structural instability. 

Fermi surfaces were calculated using a fine mesh (4800 
irreducible points, corresponding to a 48 x 48 x 32 grid- 
ding of the Brillouin zone) with a potential converged 
with 270 irreducible Brillouin zone points. Three dou- 
bly degenerate bands, which we label band 90, band 92, 
and band 94, cross the Fermi energy. The Fermi surface 
sheets formed by these three bands are shown in Figures 
|, |. In the fi gures, r, the center of the conventional 
Brillouin zone, is at the origin, i.e. at the corners of the 
cubes. The lowest energy sheet, band 90, shown in Figure 
^, consists of hole surfaces centered around F and X , giv- 
ing relatively low (^ 600 T) dHvA frequencies. The next 
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sheet, shown in Figure 0, has both hole and electron sur- 
faces stretched along the tetragonal axis. The last sheet, 
shown in Figure is predominantly an electron surface, 
although there is a large area ([001]) hole surface on that 
sheet. This sheet consists of cylinders surrounding [001] 
axis and a lattice-structure surface. The dHvA frequen- 
cies corresponding to extremal areas are given in Tables 
H and 0. 



We measured the temperature dependence of the 
dHvA amplitudes in the primary crystallographic direc- 
tions. We did not observe the anomalous temperature- 
dependence of the dHvA amplitudes noted in Ref. ||, 
where for B || [001] the amplitudes of the oscillations 
increase with increasing temperatures between 0.4 and 
1.2 K then decrease. Our measurements, carried out be- 
tween 25 and 500 mK, show decreasing amplitudes with 
increasing temperature in all cases. 

The measurements j:eported here agree substantially 
with— Cornelius et alE after a 5% correction to their 
datallj; all but one (4708 T) of the frequencies reported 
in Ref. || also are observed here. One major difference is 
the measurements along [001] where we see many more 
frequencies. Our measured frequencies in this direction 
run low as compared to calculations, but we share two 
frequencies measured by Cornelius et ai: our FV and F'g 
(see Table Also, the masses we measure for these 

orbits are in excellent agreement (see Table ^) . Cornelius 
et al. report a high mass orbit at 4686 T which we do not 
see. Another major difference is our measured masses for 
[100] are higher when compared with those of Cornelius 
et al. (see Table [V), a fact consistent with our non- 
observance of the anomalous temperature dependence. 



IV. COMPARISON OF EXPERIMENT WITH 
THEORY 

We start the comparison between the calculated fre- 
quencies and the measured frequencies with those calcu- 
lated and observed for the field along the [100] direction. 
In this case there are many more measured frequencies 
than the calculated FS gives. We attribute this situa- 
tion to strong magnetic interactions, between electrons 
on different parts of the FS, or to torque interactions due 
to the use of the cantilever. Either of these effects can 
give rise to a series of frequencies, that are combinations 
of two fundamental frequencies J13 The measured funda- 
mental frequencies in Table || are indicated as Fi through 
F7, and some of the combinations that are differences be- 
tween the fundamentals or harmonics of a fundamental 
and another fundamental frequency are shown in Figure 
|l|. It should be noted that the same fundamental fre- 
quencies are subtracted from both the fundamental 2176 
T and the first two harmonics of 2176 T (F7) to give the 
observed combination frequencies (see the inset in Fig- 
ure In one case, 1744 T is the first harmonic of the 
fundamental 874 T (Fe) frequency. 



Five frequencies ranging from 9 to 26 T obtained in 
the band calculations are too low to be observed exper- 
imentally. That is, in the field range where oscillations 
are observed (from about 7 to 18 T) there are too few 
oscillations for these frequencies to be measurable. In 
addition there are four electron orbits grouped into two 
pairs, 245 T and 215 T plus 175 T and 110 T, for which 
we observe only two frequencies, 212 T and 106 T. Again, 
the Fourier transform may not have sufficient resolution 
due to an insufficient number of oscillations to resolve all 
of these frequencies. 

The remaining five measured frequencies are in the 
same range as the predicted frequencies from the band 
calculation. Differences in detail occur, but it appears 
that for the field in this direction the measurements give 
strong support to the predicted FS. 

The agreement between the calculated and measured 
frequencies is not as good for the field in the [001] direc- 
tion. There are three calculated high frequencies, 11.268 
kT, 12.126 kT, and 12.326 kT, whereas only one fre- 
quency in this range is observed, 12.280 kT. This high fre- 
quency only is observed in the lowest temperature data, 
20 mK, and indicates a large effective mass for this fre- 
quency. The fact that only one rather than three frequen- 
cies are observed may indicate that the electron surface 
is much more cylindrical than is obtained from band the- 
ory. For the remainder of the frequencies, only the total 
number and order of the observed frequencies agree with 
the calculations. 

In the band calculation the fact that the sample is an- 
tiferromagnetic over all of the measurement field ranges 
has not been taken into account. In the model that has 
been taken for the band calculations, the Ce 4/ electron 
is completely itinerant and there are no localized mo- 
ments to cause the magnetism. Therefore, it is most 
likely that spin density waves exist in the electronic sys- 
tem and drastically alter the FS. This is the accepted sit- 
uation in CrJU and may be the case here. A much more 
involved theory is required to predict the FS under these 
circumstances. Exchange interactions could cause the FS 
to be split into separate spin up and spin down sheets. 
In order for this to be determined from dHvA measure- 
ments several hawnonics of these frequencies would have 
to be measuredJi3 We did not observe the harmonics in 
these measurements, and need to make measurements to 
higher fields at mK temperatures to do so. 



V. COMPARISON OF CERHIN5 WITH CEIRIN5 

The structure of the Fermi surface of CeRhlns is es- 
sentially the same as that calculated for Celrln5 by Haga 
et aln except that the frequency of the hole orbit per- 
pendicular to [001] arising from (what we label) band 94 
near F is 11000 T, rather than - 15000 T reported 
in Reference ^ Our own calculations for Celrlus find a 
frequency close to that calculated for CeRhlns, ^ 11600 
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The high frequency reported in Reference M was not 
seen experimentally. The difference in the two calcula- 
tions may be in our inclusion of the spin orbit interaction 
or in the fineness of the grid used to calculate the Fermi 
surface. _ 

The dHvA measurements on CelrlusEI find eight 
branches for rotations in the [100] and [110] planes of 
the tetragonal structure. Most of these branches are as- 
sociated with large quasi-2D undulating cylinders that 
show the expected l/cos(0) dependence of the measured 
frequencies with 9 being the angle at which the field is 
applied away from the [100] axis. We see a similar be- 
havior in the high frequency orbits for CeRhlns. Band 
structure calculations predict, in Ref. || and in Ref. |l^, in 
addition, that several small pieces of FS should exist in 
both CeRhlns and Celrlns. These frequencies were not 
observed in Ref. I for Celring but are observed here in 
CeRhlns . These differences likely are a result of the mea- 
surement techniques employed and not due to differences 
in the actual Fermi surfaces of the respective materials. 
Overall, the result of these angular dependent dHvA mea- 
surements on CeRhlns along with the band calculations 
show that the Fermi surface has both small pockets of 
3D character and large undulating cylinders that should 
give rise to a 2D character in other properties. 

From the temperature dependence of the signal am- 
plitudes we have calculated the effective masses of the 
observed frequencies in cases where signals were observ- 
able over the entire temperature range of measurement. 
There are two cases where no masses are reported here: 
(1) for F'l - F'3 in Table ^ the masses are sufficiently 
light that we observe no amplitude change within mea- 
surement uncertainty up to 0.5 K, and (2) for F7 in Table 
H and F'5, F'e, F'g and F'lo in Table plj the signals are 
observed only at the lowest temperatures. The masses 
are obtained from fits to the Dingle reduction factor in 
Eq. I: 



Rt ~ ttX/ sinh^TrX), 



(3) 



where A — 2TTpkT/f3B, (3 — eh/{m/'mQ)c^ and m/m^ = 
m* is the effective mass. All of the effective masses we 
obtain in this manner are greater than 1, and compare 
favorably to the values found in Ref. |^ from the high tem- 
perature dHvA data. We note that none of the measured 
masses are sufficiently large to account for the measured 
value of gamma of 400 mJ/molG=K from high temper- 
ature specific heat measurementsEl. On the other hand 
there is a loss of entropy due to magnetic ordering, and 
the inferred value of gamma below T^r is 56 mJ/mole-K, 
much more in line with the masses reported here. The 
effective masses measured on CelrlnsEi are larger than 
we observe in CeRhlus consistent with the larger value 
of gamma measured in CeIrIn5B (in which no magnetic 
order is observed). 



VI. CONCLUSION 

We have determined the electronic structure of 
CeRhlns. This determination includes both the results of 
dHvA measurements and energy band calculations. The 
major findings are: 

• The Fermi surface consists of structures arising 
from two hole bands and one electron band. The 
band calculations assume a completely itinerant 4/ 
electron from the Ce atoms and the measurements 
confirm this assumption. One of the hole surfaces 
and the electron surface are open along the c - axis 
and could give rise to some rather two-dimensional 
character to the electronic properties. 

• While the energy band calculations are for param- 
agnetic CeRhlus and the measurements are on anti- 
ferromagnetic CeRhlus , the agreement between ex- 
periment and theory is good for electron areas per- 
pendicular to the a-b plane, but deviate in size for 
orbits in this plane. This disagreement is thought 
to be due to spin density waves in the c axis di- 
rection that give rise to the anti-ferromagnetism in 
the a-b plane. 

• The electronic effective masses that are measured 
all are greater than the free electron mass and are 
anisotropic, varying from less than mo to 6mo de- 
pending on direction. 

Overall the Fermi surface of anti-fer|ipmagnetic 
CeRhlns is similar to that found previouslytl in super- 
conducting Celrln5. Thus, while their ambient pres- 
sure ground states are completely different, both of these 
materials have large electron - electron interactions and 
rather similar electronic structures. 

This work was supported in part by the National Sci- 
ence Foundation under Grant No. DMR-9971348 (Z. F.). 
A portion of this work was performed at the National 
High Magnetic Field Laboratory, which is supported by 
NSF Cooperative Agreement No. DMR-9527035 and by 
the State of Florida. Work at Los Alamos was performed 
under the auspices of the U.S. Dept. of Energy. 

TABLE I. Structural parameters of primitive tetragonal 
CeRhlns. Calculated parameters were obtained using GGA 
and brillouin zone convergence as discussed in the text. The 
parameter x is the Wykoff coordinate of displacement of 
face-centered In atoms along z. Experimental values from 
Hegger et al. 

Space Group P4/mmm (123) 



V/Ve 



1.023 



(Ve 



163.4 



c/a 

X 

B (Vea;p) 
B (VeaO 

Dos(Ef) 



1.626 
.302 
87 GPa 
78 GPa 
10.8 eV-i 



(c/aea;p 1.620) 
Xexp .306 
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TABLE II. Measured and calculated dHvA frequencies for 
CeRhlns with B along [100]. Frequencies listed for Cornelius 
et al. have been corrected by 4.5 percent (private communi- 
cation). 



calc. F (T) 



band 



F (T) 



Symb. 



F (T)El 



2139 


92 


(h) 


2176 


Ft 




849 


92 


(e) 


874 


Fe 


861 


761 


92 


(e) 


722 


Fs 


714 


551 


92 


(h) 


446 


F4 


492 


362 


90 


(h) 








342 


90 


(h) 








337 


90 


(h) 


324 


Fa 


295 


283 


94 


(e) 








245 


90 


(h) 








214 


94 


(e) 


212 


F2 


219 


175 


94 


(e) 








110 


94 


(e) 


106 


Fi 


105 


26 


92 


(h) 








22 


90 


(h) 








21 


92 


(h) 








10 


90 


(h) 








9 


90 


(h) 









TABLE V. Measured dHvA frequencies for CeRhlus with 
B along [001]. Frequencies marked with asterisks appear only 
at the lowest temperatures. Calculation of masses for the 
frequencies marked with daggers was not possible because the 
amplitudes did not change significantly over the measured 
temperature range. The frequency in parentheses is weak due 
to a low number of oscillations in the sweep range. 



Symb. 


F (T) 


m* (me) 


F (T)! 


m* (me) Ref.E 


F'lo 


12280* 








F'9 


7200* 








F'8 


6120 


6.1 ± 0.3 


6256 


6.5 ± 0.8 


FV 


3600 


4.6 ± 1.0 


3605 


4.8 ± 0.4 


F'6 


2740* 








F's 


2160* 








F'4 


1412 


4.5 ± 0.3 






F's 


562t 








F'2 


216t 








F'l 


(115)t 









TABLE III. Measured and calculated dHvA frequencies for 
CeRhlns with B along [001]. 



calc. F (T) 


band 


meas. F (T) 


Symb 


12326 


92 (e) 


12280 


F'lo 


12126 


92 (e) 






11268 


94 (h) 






6878 


92 (e) 


7200 


F'9 


5562 


94 (e) 


6120 


F'8 


4502 


94 (e) 


3600 


F't 


4331 


94 (e) 


2740 


F'6 


3951 


94 (e) 


2160 


F's 


686 


90 (h) 


1412 


F'4 


597 


90 (h) 


562 


F'3 


585 


90 (h) 


216 


F'2 


451 


90 (h) 


(115) 


F'l 



F (T) 



6256 
3605 



TABLE IV. Measured masses for dHvA frequencies for 
CeRhlns with B along [100]. Our data also are compared 



to that of Cornelius et al. 






Symb. 


F (T) 


m* (me) 


F (T) 




m* (me) Ref.E 


Fe 


874 


3.3 ± 0.6 


861 


1.31 ± 0.22 


Fb 


722 


2.0 ± 0.8 


714 


1.17 ± 0.09 


F4 


446 


3.9 ± 0.9 


492 


0.72 ± 0.14 


F3 


324 


5.5 ± 1.0 


295 


0.93 ± 0.15 


Fa 


212 


2.9 ± 0.4 


219 


0.99 ± 0.10 


Fi 


106 


3.3 ± 0.6 


105 


1.36 ± 0.23 
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FIG. 1. DFT of dHvA data taken with the apphed field 
parallel to the [100]. The inset shows some of the high fre- 
quencies that are combinations of the fundamentals. 
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FIG. 2. FFT of dHvA data taken with the applied field 
parallel to the [001]. 
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FIG. 3. Frequency as a function of angle. Calculated fre- 
quencies for [001] and [100] are indicated with bowties. 
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FIG. 4. The GGA band structure of paramagnetic 
CeRhlns calculated using experimental structural parame- 
ters. 
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FIG. 5. The GGA density of states of paramagnetic 
CeRhlns calculated using experimental structural parame- 
ters. Ce 4/ and 5d projections are shown in (a), and Rh 
and In projections in (b). 




FIG. 6. The first sheet of the paramagnetic CeRhlus 
Fermi surface. The electron surface is illuminated; these are 
hole surfaces. 




FIG. 7. The second sheet of the paramagnetic CeRhlns 
Fermi surface. The electron surface is illuminated. 
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FIG. 8. The third sheet of the paramagnetic CeRhlns 
Fermi surface. The electron surface is illuminated. 
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